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ABSTRACT 



Context. The formation of chondrules in the protoplanetary nebulae causes many questions concerning the formation process, the 
source of energy for melting the rims, and the composition of the origin material. 

/A/ms. The aim of this work is to explore the heating of the chondrule in a single precursor as is typical for radiation hydrodynamical 
shock waves. We take into account the gas-particle friction for the duration of the shock transition and calculate the heat conduction 
into the chondrules. These processes are located in the protoplanetary nebulae at a region around 2.5 AU, which is considered to be 
the most likely place of chondrule formation. The present models are a first step towards computing radiative shock waves occurring 
in a particle-rich environment. 

Methods. We calculated the shock waves using one-dimensional, time-independent equations of radiation hydrodynamics (RHD) 
involving realistic gas- and dust opacities and gas-particle friction. The evolution of spherical chondrules was followed by solving the 
heat conduction equation on an adaptive grid. 

Results. The results for the shock-heating event are consistent with the cosmochemical constraints of chondrule properties. The calcu- 
lations yield a relative narrow range for density or temperature to meet the requested heating rates of R > 10"* K h"' as extracted from 
cosmochemical constraints. Molecular gas, opacities with dust, and a protoplanetary nebula with accretion are necessary requirements 
for a fast heating process. The thermal structure in the far post-shock region is not fully consistent with experimental constraints on 
chondrule formation since the models do not include additional molecular cooling processes. 

Key words, meteorites, meteors, meteoroids - interplanetary medium - protoplanetary disks - shock waves 
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1. Introduction 

Chondrules (Greek Chondros, grain) or chondren are silicate 
globules with sizes of 0.2 to 2 mm and consist mainly of olivine 
and pyroxene. They make up 20 to 90 % of the meteorite vol- 
ume, which is why a whole group of meteorites is called chon- 
drites. The mineral structure of the chondrules together with 
their chemical composition tell a lot about their formation pro- 
cess. Based on their observed consistence and structure the heat- 
ing has to occur very quickly at a rate of > 10"^Kh"', which 
strongly favours melting of the origin material. A subsequent 
rapid cooling within a few hours is also r equired to produce the 
observed mineralogical features (see, e.g. iHood & KringI (Il996l) 
for a detail discussion of these basic properties). The analysis 
of various chondrules proves that multiple heating events have 
been present to generate their final structure. Based on these in- 
vestigations the origin of the chondrules lies within regions of 
the early protoplanetary nebulae and their age is considered to 
be between 1 and 4 mil lion years after the formation of the solar 
system (e.g. lAmeliiTer al. 2002; Kita et al. 2005). 

There exists a large number of models for chondrule for- 
mation. All of them contain more or less plausible assump- 
tions. Owing to the aforementioned structural characteristics of 
the chondrules, the duration, temperature, and pressure at the 
heating event have to occur within a spatially loc alised zone o f 
the protoplanetary nebulae as summarised, e.g. bv lBossl(ll996h . 
Shock waves appear to meet most constraints and this scenario 



is currently the most p la usible mode l , as d iscussed, e.g. by 
IConnollv & Lovd (Il998h . iDesch et al.1 (l2005l) . Since the mod- 
els based on heating in shock waves make fewer assumptions 
on the physical environment, we will focus on these events, 
which have been investigated b y a large amount of detailed 
model li ng, e.g . Hood & Horanvi ( ;199lL 1993), Ruzmaikina & Ig 
("19941), iHood (1998), Desch & Connolly (2002), Miur a et al. 
(2002), and Miura & Nakamoto (2005, 2006). Ciesla & Hood 
(2002) and ICiesla et al.l (1200 4) developed gas-particle suspen- 
sion models with significant heating via radi ation from the other 
chondrule-sized particles. In the models of iDesch & Connolly! 
(|200A the heating process is caused either by frictional heat- 
ing between gas and dust paiticles when traversing an adia- 
batic shock wave or by thermal conduction from the shock- 
heated gas to the cooler dust particles. Desch and Connolly 
also took significant heating by radiation from other chondrules 
into account. Nevertheless, without detailed modelling of the 
physical processes within the protoplanetary disc the origin of 
these shock waves remains unanswered. In the literature several 
mechanisms have been suggested to explain the shock waves, 
i.e. accretion shocks on the surface of the protoplanetary disc 
("Ruzmaikina & Ip' 1994^, infall of clumps of gas onto the neb- 
ulae (Tanakaetal. 1998|), bow s hocks produced by planetesi- 
mals on eccentric orbits (Weidenschilling et arill998l) . and X- 
ray flares from the young Sun CNakamoto et al.ll2005l) . Clearly, 
all these scenarios give only a simplified picture of the numerous 
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interactions between radiation field, and gas, and dust particles. 
In particular the treatment of adiabatic shock waves is mostly 
made with a simple equation of state and based on the equations 
of hydrodynamics without calculating a consistent radiative pre- 
cursor region. 

The recent work of iMorris & DeschI (1201 Ol) is based on a 
hydrodynamical shock model including a complete molecular 
line cooling due to H2O, a treatment of the radiation field, and 
improved opacities of the solids. These authors found the ef- 
fects of molecular cooling to be minimal because the combi- 
nation of high column densities of water, hydrogen recombina- 
tion/dissociation and radiation from the upstream regions reduce 
the rapid cooling times of chondrules in the post-shock region. 

In this work we investigate the reaction of dust particles 
to the structure of radiative shock waves. Based on full one- 
dimensional time-in dependent equations of radiation hydrody- 
namics (RHD) (e.g. iMihalas & Mihal"^ [r984f). we follow the 
evolution of a spherical dust particle that is swept through these 
shock transitions, which are parameterised by the Mach number 
In this first study we include the effects of gas-particle drift but 
the chondrules are still assumed to preserve their spherical shape 
and are not directly coupled through their radiative properties to 
the equations of RHD. Their temperature structure is obtained by 
solving a heat conduction equation allowing a thermal expansion 
of the dust particles. Thanks to the full system of RHD, the over- 
all structure of the radiative shock waves is calculated correctly 
together with the radiative precursor that heats the dust particles 
before they enter the shock jump. All models presented are lo- 
cated in a protoplanetaiy nebula at a typical distance of 2 5 AU, 
the most likely origin of chondrules (iHood & Kring|ll99a) . 

The next section deals with the basic equations describ- 
ing the interaction of dust particles, gas, and radiation in 
one-dimensional plane geometry together with the numerical 
method. In sect. 3 we present the results of several model calcu- 
lations followed by a discussion of the possible formation pro- 
cess of chondrules in radiating shock waves. 



and the equation of motion with gas pressure P and radiative 
(Eddington) flux H 
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where a:r is the so-called Rosseland-mean of the opacity, c is the 
light speed. The third equation is the gas energy equation written 
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containing the adiabatic index y, the Planck-mean opacity A-p, the 
Eddington radiation energy density J and the source function S . 
The next two equations are needed to describe the radiation 
field by writing down the Eddington moments of the radiation 
transport equation. Starting with the 0''' -moment we obtain the 
equation of the radiation energy density J 
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including the Eddington radiation pressure K. The P'-moment 
leads to the equation for the radiation flux density H, 
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This system of ODEs has to be solved for radiative shock 
waves together with additional closure conditions, and boundary 
conditions. Hence, we need an equation of state and we adopted 
in this first study an ideal gas 



A* 
The source function S is given by 



(6) 
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2. Basic equations 

2.1. Equations of radiation liydrodynamics (RHD) 

The equations of radiation hydrodynamic (RHD) treat the in- 
terplay of matter and radiation by including the momentum 
and energy exchange between the two components. A detailed 
derivation of these equations can be fo und in Landau & Li fschitz 
( Il981h and IMihalas & Mihalal ( 11984 . To investigate the physi- 
cal behaviour of dust particles we restricted the problem to time- 
independent planar shock waves in a protoplanetary nebula in 
this exploratory study. We can neglect self-gravity of the gas, and 
the resulting system of ODEs can be solved without the need of 
artificial viscosity to broaden shock waves. We furthermore as- 
sumed the Eddington-approximation as closure condition for the 
radiation field as well as an ideal molecular gas with a constant 
adiabatic index of y = 7/5 for the equation of state. 

The following three equations are the stationary, plane- 
parallel Euler equations for the gas interacting with the radiation 
field where p denotes the gas density and u the gas velocity. We 
begin with the equation of continuity 



— (Pm) = , 
dx 



(1) 



for LTE, which we assumed throughout. T is the temperature 
and cr the Stefan-Boltzmann constant. The corresponding radia- 
tion temperature Tiad can be obtained from the radiation energy 
density by 



J 



cr 



rad 



(8) 



Since the above equations contain three moments, J, K, and H of 
the grey specific radiation intensity, we need an additional clo- 
sure condition for the radiation field. Within the limit of isotropic 
intensity distributions we can use the well-known Eddington ap- 
proximation 



■' J 3 



(9) 



relating the second to the zeroth moment. For the radiating shock 
waves in a dusty environment such as the solar nebula with its 
high optical depths this approximation is usually valid. 

2.2. Boundary conditions 

According to the previous section our five basic variables are 
given by the gas density p, the gas velocity u, the gas pres- 
sure P, the radiation energy density J, and the radiation flux H. 
Hence, we have to specify five far up-stream values at x ^ -00. 
Denoting this gas density by po and the gas velocity by uq, the 
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gas pressure Pq also fixes the gas temperature Tq through the 
equation of state (|6]l. From the equation of motion (|2]l we see 
that only Hq - Q leads to consistent boundary conditions, and 
the equation of radiation energy (|4]i requires J = S and together 
with the source functions we have To = Trad.o, i-e the radiation 
temperature is equal to the gas temperature for the far up-stream 
region. The same conditions hold for the far down-stream re- 
gion jc — > +00 for stationary problems (e.g. Zel'dovich & Raizer 
Il967h . As already discussed in the literature, these conditions 
appear to be not fulfilled in all shock computations used to study 
the formation of chondrules. Because the radiative flux H is zero 
at both boundaries, a maximum (or minimum) has to occur in be- 
tween and any physical possible solution has a maximum located 
at the shock position x = 0. Ignoring for this discussion the dif- 
ference between the Rosseland- and Planck mean, we introduce 
the optical depth t by 



dr — Kp dx 



(10) 



and can transform the system of stationary RHD equations into 
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We see that the optical depth defines through (Kpy^ a length 
scale and hence (/cpowo)"' yields a typical time scale available 
for heating of particles in the precursor region. 

The properties of the shock waves are parameterised by the 
Mach number M in the up-stream region; M is given for an adi- 
abatic sound velocity ao by 



M - — with ai - y — - — To . 
ao Po H 
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A detail ed analytical treatment of RHD shock waves can be 
found in lZel'dovich & Raized (1 19671) . We used the inverse com- 
pression ratio 77 — po/p, and a value of 
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(17) 



(y + 1)M2 

has to be reached in the far down-stream region behind the shock 
wave. 

A close view on the aforementioned stationary RHD- 
equations reveals that the model assumes a single shock wave 
with a precursor for the melting of the chondrules. However, 
the complexity observed in meteorites favours formation regions 
where many precursors may exist, e.g. due to the interaction 
with bow shocks of moving planetesimals and to accretion pro- 
cesses. A more definite answer to the evolution of chondrules 
therefore requires a more elaborate treatment of the variou s in- 
teraction processes between the gas, dust particles and radiation 
fields within time-dependent flows. In particular, in another step 
it will be necessary to consider also the momentum and energy 
exchange between gas and particles because this study only takes 
into account the heating of the dust particles embedded in the 
gas without additional energy and momentum exchange. Hence, 
the cooling rates in the post shock regi on can disagree with data 
coming from furnace experiments (e.g. lHewins et al.ll2005l) . 



2.3. Gas and dust opacities 

As mentioned before, we neglected the backreaction of the dust 
particles on the gas, e.g. the evolution of the dust particle's size 
distribution function during the shock transition, which also in- 
fluences the transparency of the gas-dust mixture. The opaci- 
ties adopted in our computations were computed for evolved 
planetary discs with a chemical composition like our solar sys- 
tem, i.e. also a solid-to-gas ration of about 0.01. However, the 
chondrules may be fo rmed in a dust-enriched environment (e.g. 
ICiesla & Hoodl2002l) where this solid-to-gas ratios can be 100 or 
even 1000 times the solar ratio. In these zones the particle den- 
sities are much higher than the adopted values and the passage 
of a shock wave through this particle-gas suspensions requires 
a more detailed treatment of the opacities ( Desch et al. 20051). 
Increasing the dust-to-gas ratio increases the total opacity, and 
therefore all length scales defined through dT = Kpdx will be- 
come smaller, in particular the time tj, a chondrule needs to pass 
the radiative precursors (cf. Fig. |2) is shortened. 

Since we use frequency-integrated moments of the radiation 
field, the frequency-integrated gas and dust opacities that enter 
RHD equations, in particular the frequency-dependent opacity 
Ky, have to be integrated to obtain either the Planck-mean /cp for 
the optically thin case or the Rosseland-mean /cr in the optical 
thick case. This use of opacity tables for a dust-gas mixture can 
only be a first approximation to a more realistic situation where 
the dust particles are processed by the interaction with larger 
bodies. 

In particular, we used the Rosseland - and Planck-mean dus t 
opacities for protoplanetary discs from ISemenov et al.l (1200 3h . 
which are based on optical constants for spherical dust with 
porous aggregate particles and normal silicates. 



2.4. Gas-particle friction 

When the gas passes the shock wave, the velocity changes in 
the precursor region, the adiabatic shock, and the subsequent 
postshock cooling region. A particle changes its velocity through 
collisions with the moving gas and will therefore need some time 
to relaxate to zero relative velocity. The speed of a chondrule uc 
will be higher compared to the gas since it maintains its initial 
up-stream velocities in the absence of friction. Consequently, the 
exposition of dust grains in the hot shock zone is shorter com- 
pared to the gas. To account for this effect, we solved the fol- 
lowing system of equations for a moving chondrule at a position 
xc 

dxc dt 

Uc - or uq 1=0. (18) 

dt dxc 

The equation of motion for a chondrule can be written as 

duc dt duc 

Fc^mc—r- or Fq- otc-^— =0. (19) 

dt dxc dxc 

The simple formulae for the frictional force Fc is related to the 
relative or drift- velocity uo by 

Fc nrfjpCoUQ — 0, with ud — uc — u , (20) 

where mc is the mass of the particle and re its radius. Co is the 
dra g coefficient given for the numerical calculations according 
to (.Sandin & Hofneill2003h 



Cd - Co,difF + Cd r 



(21) 
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Fig. 1. Gas and particle velocities as a function of the distance 
from the shock wave in [km] for a Mach number M = 4.1. 
Top: Gas velocity and particle velocity (dotted line). Bottom: 
The drag force on a particle with radius re = 1 mm. 



Both terms are calculated as follows 
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s defines the type of collision, i.e. diffusive collisions lead to 
e = and specular collisions to e = 1, Tc denotes the tem- 
perature of the particles, T is the temperature of the gas and 
v„p is the thermal speed of the Maxwellian velocity distribution. 
Figure [T] shows the velocity difference between gas and particles 
and also the force on the particles plotted for a shock wave with 
Mach number M - 4.1, a preshock temperature To - 300 K, 
a preshock density po - 1.5 ■ lO^^kgm"-' and an adiabatic in- 
dex 7 = 7/5. The radius of the particle is set to re = 1 mm 
and the density corresponds to 3.2- lO^kg m"^ assuming forsterit 
Mg2Si04. 

The stopping length /, of a typical chondrule with radius re 
and density pc can be estimated by equating the mass of a dust 
particle with the gas mass contained in a cylinder of length /, and 



a cross section corresponding to the dust particle. This simple 
approximation leads roughly to 



, PC 

k - re — 

P 



(25) 



and for the above values (also used in Fig. [TJ we derived l^ ^ 
2000 km. Hence, we expect relative velocities between gas and 
chondrules of about ±2000 km in the vicinity of the shock waves 
in agreement with Fig.[T] 

2.5. Heat conduction in particles 

As a first step we studied the heating of a spherical dust particle 
by heat conduction where the temperature changes because of 
the passage through the shock at the particle surface. In the sim- 
plest approximation we assumed a spherical shape of the particle 
as an initial condition in the upstream region. The particle is also 
assumed to be in thermal equilibrium with the surrounding gas, 
to be located many optical depths away from the radiating shock 
wave, and to move with the gas, i.e. no drift velocity ud = 0. 

For calculating the particle temperature we solved the one- 
dimensional conduction equation in spherical geometry. Since 
the gas moves at high velocities and the diameter of the chon- 
drule is small (at maximum a few millimeters), we can assume 
that the temperature is homogeneous over the spherical surface. 
For further applications we rewrite the heat conduction equa- 
tion in a volume-integrated form that also allows changes of the 
radial extension or the geometrical shape for heating (and melt- 
ing). For the particle temperature Tc the spherical heat equation 
reads 



dt r^ dr \ dr 



(26) 



where D specifies the heat conduction coefficient. This coeffi- 
cient can be written as 



£) = 



PcCs 



(27) 



A is the heat conduction efficiency, pc is the chondrule density 
and Cs its specific heat. 

For the numerical solution o f this PDE we us ed a standard 
implicit difference scheme (e.g. l Richtmver '1967'). which leads 
to a tridiagonal system of equations for the temperature distribu- 
tion Tc(r) at the different particle radii. In this formulation the 
total radius can change in time due to a thermal expansion of the 
particle. The initial condition is specified by a constant temper- 
ature within the particle T(r, 0) = 7'o(r) for < r < re, and later 
on we assumed an outer boundary condition of r(rc, t) - T{xc) 
as well as dTc/dr = at the centre. The particle position xc is 
determined by integrating Eq. (fTST l. 

2.6. Numerical calculations 

For the numerical solution of the RHD system (Eqs. [T]|7]i to- 
gether with the opacity tables we adopted a standard package 
SLGA2 (Raith, Schoenauer & Glotz 1984) for ODEs, written in 
FORTRAN. We sought solutions where Hq = but dHo/dr > 0. 
Together with the remaining up-stream values the equations 
were integrated towards the shock position and matched by the 
usual Rankine-Hugoniot conditions to the down-stream values. 
A continuous transition of the gas temperature T is not possible 
but the radiation energy density J and the radiation temperature 
Trad cross the shock wave without discontinuity. 
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On top of the RHD shock solutions we followed the path of 
a dust particle by integrating Eqs. [1811241 Numerical experiments 
show that the type of gas-particle collision is negligible for our 
calculations and we adopted s - 0.5 throughout. During the tra- 
jectory the particle temperature Tc is not yet known, but should 
be in between radiation temperature T^ad and gas temperature T, 
and for these computations we considered the gas temperature 
as the appropriate outer boundary condition for the chondrules, 
i.e. Tc(x, r^rc) = T(xc). 

The calculation of the heat conduction in the dust particles 
was made by numerically solving the tridiagonal equation sys- 
tem, which we obtained from discretisation of the heat conduc- 
tion equation in spherical symmetry (u. v. 12.5b . We also included 
the thermal expansion and estimated the amount of radial in- 
crease by taking a linear expansion coefficient of a - 9.4 ■ 
10"^ K"' together with a temperature change of AT" - 1500 K. 
We obtained an increase of Arc = 1 .5 ■ 10"^ mm or about 1 .5% 
for a mm-sized chondrule. Clearly, this effect is negligible but 
we adopted a finite volume discrete version of Eq. (|26] | for sub- 
sequent applications, studying also the temporal evolution from 
fluffy to spherical particles. 

3. Results 

For all RHD shock wave models listed in Table[T]we considered 
a spherical dust particle with a radius of r = 1 mm. According 
to iHood & Krind (Il996h . this size is close to the largest chon- 
drules and we assumed a molecular gas with an adiabatic index 
of 7 = 7/5 throughout. According to Hewins et al. (2005), the 
peak temperatures during the forming process of the chondrule 
are in a range between T = 1700 - 2100 K. These gas temper- 
atures destroy all molecules and thereby also change the adia- 
batic index to a value of y = 5/3. However, to study the effects 
within the radiative precursor as well as the radiative cooling 
effects behind the shock waves, we tried to keep the equation 
of state as simple as possible for the moment. The chondrules 
were assumed to consist of forsterit (Mg2Si04) with the fol- 
lowing particle characteristics. The particle density was set to 
p = 3.2-10-'kgm"-', the heat conductivity A = 1.33 W/mK and 
the speci fic heat Cs = 1.01 kJ/kgK were taken from the litera- 
ture (e.g. lKrishnaiah. Singh & Jadhall2004tlYingweill2004l) . The 
linear heat expansion coefficient is or = 9.4- 10~* K~' . 

The protoplaneta ry disc with accretion has been modelled 
by [Bell et all ([I997h . where an angular momentum transport 
coefficient a = 10~^ leads to a mass accretion rate of M = 
1 • 10 ** Mq y"'. Outside of 5 AU these discs are gravitationally 
unstable. Taking a typical solar distance of 2.5 AU we found a 
gas density of p ~ 10"^ k g m~^ in the m id plane at a temperature 
T ~ 300 K. According to lHood & Kring ( 1996) . the temperature 
in the chondrule precursor region must not exceed T = 650 K 
so that in chondrules the existence of chemical compounds such 
as FeS can be explained. Furthermore, the melting time has to 
be in the range of about 100 s, which requires heating rates of 
R > 10"* Kh' and the subsequent te mperature decline h as to 
occur over a period of an hour or more ([Hewins et al.ltZOOSt . 

In Figure [2 a schematic curve illustrates the basic temper- 
ature structure as required from the cosmochemical evidence 
found within chondrules. Table [2] lists several models with a 
maximum temperature between Tpeak - 1700 K and Tpeak - 
2100K (see figure [3] and [4]l. The initial density po and the tem- 
perature Tq in the up-stream region {Xp, tp) correspond to dif- 
ferent Mach numbers M of the shock transitions. The adopted 
numbers vary around values suggested by simple models of the 
protoplanetary nebulae. 
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Fig. 2. Schematic temperature curve as required by chondrule 
characteristics with preheating wave To - 300 K up to T - 
1200 K in fv = 400 s {R > 10'* K h'). The shock front heats the 
gas to T+ - 1800K and during the subsequent relaxation time 
the temperature drops to Ti - 1200 K. 

Since the spatial structures within a radiative shock 
can be approximated by exponential functions, see. 



e.g. IZel'dovich & Raizen ( [19671) . we assumed for our analysis 
of the precursor region (xp, tp) and the relaxation zone (Xr, t,) 



T(x) :^ r_e 



-M/x„ 



and T{x) ^ T+e 



-M/x, 



which also defines the flow time through these regions by 



f" dx r' dx 

fp = — - and fr = — - . 

J u(x) Jo u(x) 



(28) 



(29) 



X denotes the distance from the shock wave located at jc = 0. 
The quantiles are also given in Table [T] allowing a more detailed 
characterisation of the shock structures with their typical length 
and time scales shaping the dust particles. 

The temperature difference between the surface and the cen- 
ter of a chondrule is relative small, as seen from the values of 
^T„,ax- Figure [5] shows the maximum temperature in a spherical 
body in dependence of its radius with Tpeak - 2100Kby transi- 
tion of RHD shock fronts. 



4. Discussion 

For model 1 with initial temperature Tq - 300 K and peak tem- 
perature Tpeak — 1700 K an initial density Po = 2.0-lO^^kgm"^ 
was chosen so that the requirement of a heating rate R > 
10** Kh' is fulfilled. Model 2 shows that at initial density of 
Po = 1.0- 10"'' kg m"-' the necessary heating rate is not achieved. 
A higher initial density po, like in model 3, leads to significantly 
higher heating rates. Models 4 and 5 show variations of the ini- 
tial temperature Tq, based on model 1 . A temperature decrease 
leads to higher heating rates (model 5). The necessary heating 
rates cannot be reached any more by model 4. Models 1 to 5 are 
almost subcritical shock waves, models 6 to 10 with peak tem- 
peratures Tpeak = 2100 K constitute supercritical shock waves. 
With density and temperature variations models 6 to 10 react in 
the same way as models 1 to 5, but the changes have a greater 
effect on the duration of the preheating wave and thereby on the 
heating rate. The increase of temperature in supercritical shock 
waves at the beginning of the preheating waves fulfills the nec- 
essary heating rates. The particles stay for a long time in the 
preheating wave until the shock occurs. 



Table 1. Parameter of RHD shock waves with y - 1 15 
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Model 


1 


2 


3 


4 


5 


po[lO-"'kgm-1 


2.0 


1.0 


3.0 


2.0 


2.0 


ro[K] 


300 


500 


100 


r+[K] 


1700 


M 


3.97 


3.97 


3.98 


2.92 


7.23 


Xp [km] 


83 


3600 


370 


1100 


980 


Xr [km] 


175 


225 


132 


220 


131 


t„ [s] 


179 


780 


79 


240 


199 


?,[s] 


43 


55 


32 


60 


30 


R{\{f Kh-'] 


-2.4 


-0.6 


-5.8 


-0.9 


-5.3 


Ar,„„, [K] 


13.2 


12.8 


13.4 


9.0 


17.2 




Model 


6 


7 


8 


9 


10 


Po[10-''kgm-'] 


5.83 


5.7 


6.0 


5.83 


ro[K] 


300 


400 


200 


r+[K] 


2100 


M 


4.48 


4.48 


4.48 


3.80 


5.61 


Xp [km] 


1150 


5600 


310 


76000 


106 


X, [km] 


122 


122 


122 


125 


119 


tp [s] 


232 


1220 


59 


17100 


20 


t, [s] 


36 


36 


36 


38 


32 


R[WKh-'] 


-2.2 


-0.4 


-8.1 


-0.03 


-30.7 


AT-,,,,,, [K] 


26.8 


26.2 


28.8 


20.7 


33.7 



Investigating the range of possible initial conditions as con- 
strained by models of protoplanetary discs, we found that the 
gas can be heated up to Tp^ak - 2100 K and that the Rosseland- 
and Planck-mean opacities with dust can meet the basic require- 
ments of shock-based models. These RHD shock waves provide 
the requested heating rates. Calculations using a monoatomic 
gas with y - 5/3, opacities without dust (Alexander opacities of 
iGrevesse & Noels (1993)) and a disc models without accretion 
(e.g. the minimum mass nebulae of B lack & Metthews (1985 )) 
reveal that the particles are propagating within the preheating 
wave for days to years, which does not match the required heat- 
ing rates. Molecular gas and the dust-enriched opacities result 
in better agreements with the observational facts of chondrule 
rese arch. The model s of a protoplanetary disc including accre- 
tion dBell et alJll997h have much higher gas densities in the mid 
plane around 2.5 AU. This gas densities are essential to compute 
the requested heating rates. Summarising we emphasize that the 
increase of the up-stream density po together with a decrease of 
the temperature To makes molecules and therefore y - 7/5 as 
well as higher opacities possible. All these effects reduce the du- 
ration of the dust particle in the preheating wave and increase 
the heating rates. 

The calculated temperature differences between surface and 
core of the chondrules (Fig. |5]l are correct as long as the chon- 
drules remain solid particles. But up to now we have not incorpo- 
rated phase transitions in this chondrule model. Clearly, strong 
shock waves with high Mach numbers M as well as shock waves 
with smaller preheating waves will produce larger temperature 
differences but are unlikely to ensure the stringent cosmochem- 
ical constraints. As depicted in model 10, the preheating wave 
vanishes completely and the sudden increase of the temperature 
at the outer boundary yields the largest temperature differences 
within chondrules. For these conditions it is easier to explain 
partial or shell melting of the particles. 

Immediately after the shock passage we observe a rapid 
cooling on a time scale of minutes owing to radiative losses de- 
termined by the opacities. This effect is followed by a longer 




Fig. 3. The temporal structure of the temperature within a RHD 
shock wave reaching a maximum of Tp^ak - 1700 K in units of 
hours relative to the passage through the shock wave. The top 
panel plots the variations in pQ, Tq - 300K = const; bottom: 
Variations in Tq, po - 2-10"^ kg m"^ - const; solid line: model 1, 
dotted line: model 2, dashed line: model 3, dash-dot line: model 
4, dash-dot-dot line: model 5 



time scale of several days where non-equilibrium processes fur- 
ther reduce the post shock temperature of the gas and the chon- 
drules. The rapid cooling rates of our models of about 10'*K/hr 
seem to disagree with furnace experiments, which indicate val- 
ues of around lOOK/hr The reason for this discrepancy is due 
to different time scales involved during cool ing processe s in the 
post shock region. The calculations of, e.g. iDesch et all (l2005h 
include the energy and momentum exchange between the chon- 
drules and the gas. Therefore the temperature structure in the 
post shock region is controlled over longer time scales by an 
additional cooling term A(x), which includes the radiative cool- 
ing by HaO-molecules. Consequently, the far down-stream gas 
temperature can be set to the far up-stream temperature Tq. This 
leads to overall cooling rates of lOOK/hr, which cannot be ob- 
tained within our physical description, which neglects this addi- 
tional cooling term. The boundary conditio n demands Tn = T i 
for the gas temperature as specified, e.g. in iDesch et al.l (l2005h . 
and the post shock region cools down to the initial tempera- 
ture by radiative looses into the surrounding protosolar nebula. 
However, this additional cooling occurring on a longer spatial 
and temporal scale reduces the overall cooling rates to less than 
lOOK/hr From furnace experiments (Hewins et al. 2005) it has 
become clear that the exact shape of cooling path has no influ- 
ence on the final coolin g rate. 

iMorris & DeschI (1201 Ob reveals how ad- 



The recent work of 



ditional physical mechanisms like molecular line cooling due 
to water influence the cooling history in the post-shock region. 
However, these authors emphasize that the post-shock region is 
only slightly changed because several effects partly cancel each 
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Fig. 4. Heating to Tp^ak = 2100 K within RHD Shock waves; top: 
variations in p^), Tq = 300K - const; bottom: variations in To, 
po = 5.83-10"''kgm"-^ = const; solid line: model 6, dotted line: 
model 7, dashed line: model 8, dash-dot line: model 9, dash-dot- 
dot line: model 10 
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Fig. 5. Maximum relative temperature difference within a spher- 
ical chondrule as a function of its radius in [mm] for shock waves 
with Tpeak = 2100 K as summarised in Table[T] solid line: model 
6, dotted line: model 7, dashed line: model 8, dash-dot line: 
model 9, dash-dot-dot line: model 10 



other out, i.e. radiative heating from the shock front as well 
as recombination and dissociation of molecules. Summarising 
these recent results of Morris & Desch (2010) we note that the 
improved nebular shock models can meet all cosmochemical 
constraints of condmle formation. The shock speeds necessary 
to melt chondrules have to be increased, e.g. from 7km/s to 
8km/s. 

Several models have been developed during the past years to 
calculate the behaviour of dust particles within shock waves (e.g. 
iDesch & ConnoUvll2002t iMiura & Nakamotdl2005l) . However, 



none of the proposed models have included the full set of RHD 
equations with the most recent gas and dust opacities, which are 
essential for the correct structure of the precursor region . This 
zone is shaped by the radiation transmitted from the shocked 
down-stream gas and the dust particles are advected through this 
region towards the shock wave. As inferred from the cosmo- 
chemical evidence, this region provides a basic heating process 
and will therefore restrict the possible Mach numbers. As seen 
from the models (cf. Table [TJ, already moderate Mach numbers 
around M ^ 4 (for T - 300 K) meet the constraints of chon- 
drule formation with parameters reasonably within protoplane- 
tary nebula models. 

Opacities changes due to dust-enriched envir onments caused 
by fragmentation and collisional debris (e.g. ICiesla & HoodI 
i2002i) may significantly increase solid-to-gas ratios, and val- 
ues of 100 or even 1000 times the solar ratio can be ex- 
pected. Clearly, shock waves propagating in such a densely dust- 
populated enviro!lSl£SlI23Hi'"^s ^ more detailed treatment of the 
opacities (iDesch et al.ll2005l) . as adopted in these simulations. In 
particular, all time- and length scales controlled by the opacity 
are then changed. 

As mentioned already, the shape and/or the mass of the par- 
ticles can be modified as they propagate through the shock re- 
gion. First, the partial or total melting could lead to a mass 
loss by evaporation or, depending of the surface tension, also 
to a rearrangement of the surface, i.e. the particles can become 
more spherical. Secondly, drift velocities acting on the surface 
can result in an erosion and thereby decrease the particles mass. 
Thirdly, the radiative interactions with the dust particles together 
with the additional heat needed to melt the particle should be in- 
cluded into these models of radiative shock-heated chondrules. 
The dissociation of molecules also demands a more realistic 
treatment of the equation of state. These structural effects of the 
dust particles will be discussed and included in a forthcoming 
paper in more detail. 
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